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Abstract 



Chiral symmetry is dynamically broken in quenched, ladder QED at weak 
gauge couplings when an external magnetic field is present. In this paper, we 
show that chiral symmetry is restored above a critical chemical potential and 
the corresponding phase transition is of first order. In contrast, the chiral 
symmetry restoration at high temperatures (and at zero chemical potential) 
is a second order phase transition. 







Do external fields affect the symmetry properties of the vacuum? |J The answer is yes; 
in particular, it was found that, in the quenched, ladder approximation of QED, chiral 
symmetry is dynamically broken at weak gauge couplings when a uniform magnetic field is 
present Subsequently, in a more detailed exposition we show that chiral symmetry 

is restored above a critical temperature which is of the order of m , the dynamical fermion 
mass generated at zero temperature and zero chemical potential (see also [f|). In this paper, 
we consider the effects of a chemical potential on the chiral symmetry breaking in QED 
induced by a magnetic field. We find that chiral symmetry is restored above a critical 
chemical potential fi c which is proportional to with the proportional factor a function of 
the gauge coupling constant. We also find that the dynamical fermion mass m M generated at 
nonzero chemical potential increases from m as the chemical potential [i is increased from 
zero toward fi c and exhibits a discontinuous jump to m M = when the critical chemical 
potential is crossed. The order parameter for the phase transition, the chiral condensate 
(i/jiI>), has the similar discontinuous behavior at fi c . This indicates that the corresponding 
phase transition is of first order. For comparison, we study the nature of the chiral phase 
transition at high temperature but with zero chemical potential. This phase transition is 
of second order since the dynamical fermion thermal mass as well as the corresponding 
chiral condensate approach zero continuously as the temperature approaches the critical 
temperature from below. 

Unless otherwise noted, all notations are the same as those in Ref. 0], hereafter referred 
to as I. 

The effects of a chemical potential /i can be incorporated into the study of chiral sym- 
metry breaking in an external magnetic field at nonzero temperature j|] by making the 
substitution |J: 

p -> mT(2l + 1) - /x, Z = 0, ±1, ±2, ... (1) 

for the fermion energy. The thermal photon propagator remains unchanged. One can 
easily show that the presence of a chemical potential does not modify the orthogonality 
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and completeness conditions |7]] of the basis eigenf unctions used in our formalism. It is 
then straightforward to follow the procedures described in I to obtain the gap equation for 
nonzero temperature and nonzero chemical potential: 
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Ql + 4tt 2 T 2 / 2 Q 2 + [ttT(2Z - 1) - i/i] 2 
where Q 2 = g 2 + m 2 (T, /i), Q 2 . = <?3 + 2|eif|<^, and m{T,n) is the infrared dynamical 
fermion mass as a function of both the temperature and the chemical potential. 

We sum over / in Eq.(^J) using the Poisson sum formula. The gap equation now reads 
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In the 11 — limit, this equation reduces to the nonzero temperature gap equation obtained 
in I §. 

Following the procedure described in Appendix B of I, one can obtain the chiral conden- 
sate at nonzero temperature and nonzero chemical potential as 
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where we have cut off the momentum integral at J\eH\ (because the dynamical fermion 



mass would be vanishingly small at momenta above J\eH\) and ignored the momentum 
dependence in m. The sum can again be evaluated by means of the Poisson sum formula to 
yield 
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In the limit T = /i = 0, Soo equals 1 and we get (m = m(0, 0)) 
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in agreement with the result given by Eq.(B4) in I. 

To isolate the effects of the chemical potential on chiral symmetry breaking in the pres- 
ence of an external magnetic field, we now consider the T = limit of the above gap 
equation. We note the following simplification in this limit: 

coth(^) -> 1, 

tanhpi±^) -> e{n) + e{-fi) [B{Qi + fi)-e{-fji-Qi)], 

tanh(^) -> 0(n) [e(Q 1 -fJt)-e{p-Q 1 )] + 0(-fM). 
Eq.([3]) is then reduced to 
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As a check, in the 11 = limit, we recover Eq.(65) in I, which is the gap equation for the 
case of T = and \i = 0. Notice that Eq.(^) as well as Eq.(|5]) do not depend on the 
sign of /1. For ease of writing, henceforth we drop the magnitude notation on \fi\ with the 
understanding that fi stands for the nonnegative \fi\. 

In the small /1 limit where fi -C mo, we may treat the chemical potential effects as a 
perturbation and write m 2 (0, /a) ~ m\ + 5m 2 with Sm 2 <C m\. This allows us to obtain 
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and Q 2 = q 2 + Note that Ji and I2 are the same integrals as defined in Eq.(67) of I. 
They are both finite and positive so that dm 2 ^ is positive. 

In the large \x limit where fi ^> m , we may rewrite Eq.(0) as 
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Since both g 3 -integrals are dominated by values of g 3 near the respective lower limit of 
integration, it is clear that the first term in the square bracket is larger than the second 
term. Consequently the right hand side of the gap equation is negative and there is no 
solution for large [l. This implies that there is a critical chemical potential above which 
chiral symmetry is restored. Indeed, this expectation is confirmed by the result of our 
numerical study discussed below. 

Being unable to solve the gap equation, Eq.(f|), analytically, we seek a numerical solution 
instead. It is convenient to introduce the following dimensionless variables: 

qs = T= t nrin = 7= ft = 7= Qi,a = 7= (12) 
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In terms of these variables, Eq.(H) becomes 



2a 

1 - — Iim^p), (13) 



7T 



I poo ^ poo _ 2 e -9x 



2 



where 

j /*oo /*oo 

I{m^p) = - dq 3 dq^ 

4 Jo Jo QiQ 

■ U 1 + H® 1 7® - 9{ll 7 ) (14) 

For various values of /t and , one can perform the double integration numerically. The 
result is shown in Fig. 1 where we have plotted I(m M ,/2) as a function of m M for five 
representative values of successively increasing fi from fi = to p, — 0.004243. The 
numerical solution for m M for each value of (1 can be read off by equating J(m M ,/t) to 



it /2a. For illustrative purpose, we have used a = tt/30 in Fig. 1. Going from Fig. lb to 
Fig. lc, as p, increases, one finds two solutions for rh^ instead of one. As is clear from Fig. 
la,b,c, the requirement of continuity to the fx = solution selects the larger value of 
as the consistent solution ||. Increasing fi to 0.004101, one finds again only one solution 
for rhfj,, as shown in Fig. Id. Further increase in p, results in no solution as indicated in 
Fig. le. We therefore conclude that the critical value of fx is 0.004101 (for a = 7r/30) 
above which chiral symmetry is restored. 

In Fig. 2, the dimensionless dynamical fermion mass m M is plotted as a function of 
the dimensionless chemical potential fx for a = tt/30. We note that the phase transition 
across the critical chemical potential is of first order as the dynamical fermion mass has a 
discontinuity at the transition point. The order parameter also exhibits such a discontinuity. 
The chiral condensate at T = and fx ^ can be obtained from Eq.flSl) by noting that 



Sofi = d(Jp3 + mft — a0, where m M denotes m(0, fx). We find 
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Since m M = above the critical chemical potential and > fx for fx < fx c (see Fig. 2), 
only the first integral on the right hand side of Eq. ([15]) contributes. The result is 



WW ^ - ^m M ln(^) (16) 



Clearly, the result in Eq.(^) is recovered when fx = 0. It is also apparent that the disconti- 
nuity in is reflected as a discontinuity in (f/^W a ^ * ne critical chemical potential. 

One can numerically compute the critical chemical potential fx c for different values of a. 
The result, as depicted by the solid curve in Fig. 3, shows that fi c increases with a. This 
is an expected behavior because it requires a larger chemical potential to destabilize a more 
strongly coupled chiral condensate. 

A closer examination of Fig. 1 shows that, for each value of /}, /(m M ,/t) peaks at 
= fx. One can understand this feature qualitatively. For a given value of fx, let us 



consider the integral J(m^,/t) as a function of rh^. For rh^ > jl, Eq.(|T4]) becomes 
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Clearly, I(m^,jl) increases monotonically as rh^ is decreased; so in the region where 
rh^ > jl, J(m M ,/i) reaches its maximum at m M = jl. For rh^ < jl, we may examine the 
slope 
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Note that it is zero at m M = and goes to infinity as m M approaches jl. Furthermore, ^ , . 
is positive at = and the positive term in the integrand (the term proportional to the 
5-function) increases with increasing rh^ while all the negative terms decrease in magnitude. 
In other words, the slope is always positive for m M < jl. Hence, for fixed jl, 7(rh M ,/t) is a 
monotonically increasing function of rh^ in the region m M < jl. The discontinuity of the 
slope at = jl indicates that J(m M ,/i) should peak sharply there, as we see from Fig. 1. 
We shall exploit this fact to estimate analytically the critical chemical potential. 

As the chemical potential reaches \x c , I{m^,jl c ) peaks at m M = fi c and I(jl c ,pi c ) = n/2a 
satisfies the gap equation, Eq. (|TB|) . This implies that 

= fie) = Ac" ( 19 ) 

Adopting the perturbative result in Eq.(^), we find 
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This estimate of ji c is shown as the dotted curve in Fig. 3 where the approximate analytical 
solution for mo (see Eq.(53) in I) has been used. The good agreement with the exact 
numerical result (solid curve in Fig. 3) indicates that the analytical estimate of the critical 
chemical potential is quite robust for small values of a. 

As noted above, the phase transition across the critical chemical potential is of first 
order. In I where we studied chiral symmetry breaking in an external magnetic field at 
nonzero temperature but with zero chemical potential, we showed that the chiral symmetry 
is restored above a critical temperature T c . In the remainder of this paper, we would like 
to examine this phase transition across T c in greater detail. With the introduction of the 
dimensionless temperature and dynamical fermion thermal mass: 
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the finite temperature gap equation can be obtained by setting fi = in Eq. 
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We have solved the above equation numerically. The result is shown in Fig. 4 where, for 
illustrative purpose, we have again used a = n/30. The critical temperature at which the 
dynamical fermion mass vanishes is T c pa 0.00205 for this value of a. As shown in Fig. 4, 
starting at T — T c , rfiT increases as T decreases and it remains zero as T increases indicating 
that the corresponding phase transition is of second order. 
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To find out how the dynamical fermion mass varies with temperature near (and below) 
T c , we expand J(rhr,T) around m\ = and T = T c so that for small T — T c and small 
TTi j^, Eq.(p2|) takes the form 



2a 
1 ~ — 

7T 



(24) 



It can be shown both analytically and numerically that 9J ^?^ and dJ ^z c " > are finite and 
negative. We can understand qualitatively why these two partial derivatives must have 
the same sign. According to the result in Fig. 4, increasing the temperature reduces the 
dynamical fermion thermal mass. It follows that raising the temperature must have the 
opposite effect on J(rfiT,T) as lowering the fermion thermal mass in order that the gap 
equation, Eq. (|22|) , remains satisfied. Since rh 2 -. = for T = T c solves the gap equation, we 
obtain from Eq.fl2"4"D that 

<9J(0,T C U / f\ dJ(0,f c ) . 

= K -4-^T c 1 - — + V ' J m 2 T) 25 

&T V TJ dm 2 T v ' 

which indicates that txit has the behavior 

m T ~ (l - |) 2 (26) 

as T approaches T c from below. This translates into a similar behavior for the order param- 
eter of this phase transition. 

The chiral condensate at nonzero temperature and zero chemical potential can be gleaned 
from Eq.(||) to be 



\eH\ rVW dp 3 ( ^Vz + m 2 T 
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where ttit denotes m(T, 0) and the last result is obtained by noting that the integral is 
dominated by large values of p 3 for which the hyperbolic tangent function in the integrand 
is approximately 1 (for T < T c ). Eq.(p6|) and Eq.(p7|) show that 
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as T — > T~ . 

To summarize, we have considered the separate effects of a chemical potential and tem- 
perature on the chiral symmetry breaking induced by an external uniform magnetic field in 
quenched, ladder QED. We have done so in order to highlight the effects of each. Of course, 
it is useful and interesting to consider their combined effects. Work along this direction is 
in progress. 
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FIGURE CAPTIONS 



Fig. la: J(m^,/t) vs. m M for fi = 0. J(rh^,/t) peaks at m M = /t = 0. The solution to 
the gap equation is at = 0.003790 for a = 7r/30. 

Fig. lb: J(m M ,/2) vs. m M for /} = 0.001414. Iirh^pi) peaks at = /t = 0.001414. 
The solution to the gap equation is at m M = 0.003817 for a = 7r/30. 

Fig. lc: J(m^,/t) vs. m M for \x = 0.002828. /(m^/i) peaks at = fx = 0.002828. 
Note that there are two possible solutions to the gap equation at m M = 0.002715 and 
0.003889 for a = vr/30. 

Fig. Id: I(m^,£i) vs. m M for ft = 0.004101. I(m^,(i) peaks at = fi = 0.004101. 
There is one solution to the gap equation at m M = fi = 0.004101 for a = 7r/30. 

Fig. le: I(m M ,/i) vs. rh^ for fi = 0.004243. I(m^,jl) peaks at = fi = 0.004243. 
There is no solution to the gap equation for a = 7r/30. 

Fig. 2: The dynamical fermion mass as a function of the chemical potential fi 



(both measured in units of J2\eH\) for a = 7r/30. The corresponding dimensionless critical 
chemical potential fi c is equal to 0.004101. 



Fig. 3: The critical chemical potential (measured in units of y2\eH\ ) as a function of the 
gauge coupling a. Solid curve: from exact numerical calculations. Dotted curve: from the 
estimate in Eq.(|20]). 



Fig. 4: The dynamical fermion thermal mass tjit as a function of the temperature (both 



measured in units of y2|eif|) for a = 7r/30. 
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Fig. 1b 
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Fig. 1c 
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Fig. 1e 
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Fig. 3 
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